Interaction-based quantum metrology showing scaling 
beyond the Heisenberg limit. 

M. Napolitano\ M. Koschorreck\ B. Dubost^ ^ N. Behbood\ R. J. Sewell^ & M. W. Mitchell^ 

^ICFO-Institut de Ciencies Fotoniques, Mediterranean Technology Park, 08860 Castelldefels 
(Barcelona), Spain. 

^Laboratoire Materiaux et Phenomenes Quantiques, Universite Paris Diderot et CNRS, 
UMR 7162, Bat Condorcet, 75205 Paris Cedex 13, France. 

Quantum metrology studies the use of entanglement and other quantum resources to im- 
prove precision measurement^. An interferometer using independent particles to mea- 
sure a parameter X can achieve at best the "standard quantum limit" (SQL) of sensitivity 
SX oc A^"^/^. The same interferometer^ using entangled particles can achieve in princi- 
ple the "Heisenberg limit" 5X oc A^"^, using exotic states^. Recent theoretical work argues 
that interactions among particles may be a valuable resource for quantum metrology, allow- 
ing scaling beyond the Heisenberg limiP*^. Specifically, a fc-particle interaction will produce 
sensitivity 5X oc A^"^ with appropriate entangled states and 5X oc A^-(^~i/2) even without 
entanglement. Here we demonstrate this "super-Heisenberg" scaling in a nonlinear, non- 
destructive^ ^ measurement of the magnetisation^^ of an atomic ensemble^^. We use fast 
optical nonlinearities to generate a pairwise photon-photon interactiorP^ (k = 2) while pre- 
serving quantum-noise-limited performanccEHMj, to produce 5X oc A^"^/^. We observe super- 
Heisenberg scaling over two orders of magnitude in A^, limited at large A^ by higher-order 
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nonlinear effects, in good agreement with theorjff^. For a measurement of limited duration, 
super-Heisenberg scaling allows the nonlinear measurement to overtake in sensitivity a com- 
parable linear measurement with the same number of photons. In other scenarios, however, 
higher-order nonlinearities prevent this crossover from occurring, reflecting the subtle rela- 
tionship of scaling to sensitivity in nonlinear systems. This work shows that inter-particle 
interactions can improve sensitivity in a quantum-limited measurement, and introduces a 
fundamentally new resource for quantum metrology. 

The best instruments are interferometric in nature, and operate according to the laws of 
quantum mechanics. A collection of particles, e.g., photons or atoms, is prepared in a superposition 
state, allowed to evolve under the action of a Hamiltonian containing an unknown parameter X, 
and measured in agreement with quantum measurement theory. The complementarity of quantum 
measurement^ determines the ultimate sensitivity of these instruments. 

Here we describe polarisation interferometry, used for example in optical magnetometry to 
detect atomic magnetisatioiinHHiZl similar theory describes other interferometers^. A collection of 
photons, with circular plus/minus polarisations |+) , |— ) is described by single-photon Stokes 
operators Si = |(|+),|— |)^, where the gi are the Pauli matrices and ctq is the identity. 
In traditional quantum metrology, a Hamiltonian of the form H = hX XljLi ^^^^ uniformly and 
independently couples the photons to X, the parameter to be measurecP. If the input state consists 
of independent photons, the possible precision scales as (^A' oc A"^/^, the shot-noise or standard 
quantum limit (SQL). The N~^l^ factor reflects the statistical averaging of independent results. 
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In contrast, entangled states can be highly, even perfectly, correlated, giving precision limited by 
5X oc A^"\ the Heisenberg limit (HL). 

The above Hamiltonian is conveniently written H = hXSz, where Si = J2j=i ^i^"^ ^ 
lective variable describing the net polarisation of the photons. The independence of the photons 
manifests itself in the linearity of this Hamiltonian. Recently Boixo et al. have shown that inter- 
actions among particles, or equivalently nonlinear Hamiltonians, can contribute to measurement 
sensitivity and give scaling beyond the Heisenberg limit"^. For example, a Hamiltonian H = hXS^, 
i.e., with a /c-order nonlinearity in S, contains /c -photon interaction terms s^J^^ ® si^^^ ... s^z^\ 
The number of such terms, and thus the signal strength, grows as A^^, while the quantum noise from 
the input states is unchanged. As a result, a sensitivity limit of ^A' oc A^"^ applies when entangle- 
ment is used, and 5X oc A^~(^~i/2) in the absence of entanglement'^. For > 2, this already gives a 
scaling better than the Heisenberg limit, so-called "super-Heisenberg" (SH) scaling^. Note that in- 
teractions and entanglement are compatible and both improve the scaling. The predicted advantage 
applies generally to quantum interferometry, and proposed mechanisms to produce metrologically- 
relevant interactions include Kerr nonlinearitie^, cold collisions in condensed atomic gaseP, 
Duffing nonlinearity in nano-mechanical resonator^ and a two-pass effective nonlinearity with 
an atomic ensemble^^. Topological excitations in nonlinear systems may also give advantageous 
scaling. 

In this Letter, we study interaction-based quantum metrology using unentangled probe par- 
ticles. One challenge in demonstrating SH scaling is to engineer a suitable nonlinear Hamiltonian. 
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Some nonlinearities have been shown to be intrinsically nois)P^ while others give SH scaling but 
fall short of the ideal A^~(^~i/2) under realistic condition^*^. We use a cold atomic ensemble as a 
light-matter quantum interfaceplto produce quantum-noise-limited interactions and a Hamiltonian 
of the form H = fiXS^S^) = hXSzN/2. This Hamiltonian gives a polarisation rotation growing 
with the photon number, without increasing quantum noisd^. The experiment, shown schemati- 
cally in Fig. [ij uses pulses of near-resonant light to measure the collective spin F of an ensemble 
of Na ^ 10^ cold rubidium-87 atoms, probed on the 55^1/2 5P3/2 D2 line. The experimental 
system is described in detail in the reference^SHH. The on-axis atomic magnetisation ( ), which 
plays the role of X in this measurement, is prepared in the initial state { F^) = Na by optical 
pumping with resonant circularly polarised light propagating along the trap axis z. A weak on-axis 
magnetic field is applied to preserve during the measurements. 

Pulses of Sx polarised, but not entangled, photons pass through the ensemble and experi- 
ence an optical rotation proportional to { F^). The light-atom interaction Hamiltonian ^eff = 
a^'^'^FzSz + (^^^^FzSzN /2 describes this paramagnetic Faraday rotation^^. Both the linear term 
a^^^FzSz and the nonlinear term P'^^^FzSzN /2 cause rotation of the plane of polarisation from Sx 
(vertical) toward Sy (diagonal). Detection of Sy then allows estimation of F^. As described in the 
Supplementary Information, a^^^ and (5'^^^ depend on the optical detuning A relative to the F = 
1 ^ = transition, in particular ^(^^(Ao) = for the specific detuning Aq 27r(468.5 MHz), 
allowing a purely nonlinear estimation to be studied. 

The rotation angle is = ( ) [A(A) + B{/\)N]/2 where A oc a^^^ and B oc (3^^^ account 
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for the temporal pulse shape and geometric overlap between the atomic density and the spatial 
mode of the probe. The shot-noise limited uncertainty in the rotation angle, due to quantum uncer- 
tainty in the initial angle, is Sep = A^"^/^/2. A contribution ( ) B{A)5N/2 from initial number 
fluctuations SN = (A^)"^^^ is negligible for small rotation angles. This gives a measurement 
uncertainty 

Scj) 1 

indicating a transition from SQL scaling SF^ oc N~^/^ to SH scaling SF^ oc A^~^/^ with increasing 

Two regimes of probing are used: the linear probe consists of forty 1 |lis pulses (total illu- 
mination time tl = 40|lis) spread over 400|lis with detuning Al 3> Aq. This gives A Ni^B, 
i.e., linear estimation and, as described by Koschorreck et alW, provides a projection-noise-limited 
quantum-non-demolition (QND) measurement^"^ of F^, with uncertainty at the parts-per-thousand 
level^. The nonlinear probe consists of a single tnl = 54 ns FWHM, Gaussian-shaped, high- 
intensity pulse with A^NL photons and detuning Aq, so that A <C N^lB. Crucially, having two 
probes allows us to precisely calibrate the nonlinear measurement using a highly sensitive and 
well characterised independent measurement of the same sample. 

We probe the same sample three times for each preparation: First with the linear probe, 
which gives a precise and non-destructive measurement of ( ) via a rotation Then with the 
nonlinear probe, contributing with a rotation which is calibrated against the "true" value (i.e., 
with negligible error) provided by the previous linear probe. Third, a second linear probe is used 
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to estimate the damage to the atomic magnetisation 77 = 1 — caused by the nonhnear probe. 

The hnear probe is cahbrated using quantitative absorption imaging to measure A^^, and we 
find ^(Al) = 3.3(1) X 10"^rad per atom. The cahbration of the nonhnear probe against the 
first hnear probe is shown in Fig. |2} We repeat the above pump/probe sequence while varying 
Na in the range A^^ = 1.5 x 10^ to A^^ = 3.5 x 10^ to generate a 0nl correlation plot for 
a given A^nl- Since both and (/)nl are linear in A^^, we use linear regression to find the slope 
b = d(f)^i^/d(f)i^ = 5(Ao)A^nl/^(Al) for that value of A^nl- The experiment is repeated varying 
the number of photons A^nl in the nonlinear pulse. 

The observed h vs. A^nl, shown in Fig. [2^, is well fit by a simple model including saturation 
of the nonlinear response: 



6nl ^ ^(Ao)A^NL 1 

^(A^) 1 + N^jNt'^ ' 



with a saturation parameter N^^^ = 6.0(8) x 10^ and the nonlinear coupling strength 5(Ao) = 
3.8(2) X 10"^^rad per atom per photon. 

The noise in the nonlinear probe, again as a function of A^nl? is determined from the 0nl 
correlation plots. As illustrated in Fig.|2|3-c, the residual standard deviation of the fits indicates the 
observed uncertainty A(/)nl, which includes the intrinsic uncertainty 5(/)nl and a small contribution 
from electronic noise. In Fig. [s we plot the fractional sensitivity 5Fz^^^ / { ) vs. A^nl, calculated 
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using equation ^ and considering the whole polarised ensemble, ( ) = 7 x 10^. In agreement 
with equation ([l]), the log-log slope indicates the scaling 5Fz^^'' oc A^^l^^ within experimental 
uncertainties in the range A^nl = 10^ to A^nl = 10^, and SH scaling, i.e., steeper than A^"^, over 
two orders of magnitude A^nl = 5 x 10^ to A^nl = 5 x 10^. 

Results of numerical modelling using the Maxwell-Bloch equations to describe the non- 
linear light-atom interaction are also shown in Fig. |3j Two curves are shown, for detunings 
Aq ± 27r(200 kHz), covering the combined uncertainty in A due to the probe laser linewidth and 
inhomogeneous light shifts in the optical dipole trap. As expected from equation ([T]), this alters the 
sensitivity only at low A^nl- The model is described in detail in the Supplementary Information. 

For photon numbers above A^nl ^ 2 x 10^, the saturation of the nonlinear rotation alters the 
slope. This can be understood as optical pumping of atoms into states other than \F = I.ttif = 1) 
by the nonlinear probe. The damage to the atomic magnetisation = 1 — (/>lV^l, shown in Fig. [5] 
remains small, confirming the non-destructive nature of the measurement. The finite damage even 
for small A^nl is possibly due to stray light and/or magnetic fields disturbing the atoms during 
the 20 ms period between the two linear measurements. At large A^, high-order nonlinear effects 
including optical pumping limit the range of SH scaling. 

The experimental results illustrate the subtle relationship of scaling to sensitivity in a non- 
linear system. For an ideal nonlinear measurement, the improved scaling would guarantee better 
absolute sensitivity for sufficiently large A^. Indeed, when the measurement bandwidth is taken 
into account, the nonlinear probe overtakes the linear one at A^ = 3.2 x 10^ where both achieve a 
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sensitivity of 1.1 x 10^ spins Hz"^/^. As a consequence, the nonlinear technique performs better 
in fast measurements. In contrast, when measurement time is not a hmited resource, the com- 
parison can be made on a "sensitivity-per-measurement" basis, and the ideal crossover point of 
3.2 X 10^ spins at = 8.7 x 10^ is never actually reached, due to the higher-order nonlinearities. 
Evidently SH scaling enables but does not guarantee enhanced sensitivity: for the nonlinear to 
overtake the linear, it is also necessary that the scaling extend to large enough A^. The comparison 
shows also that resource constraints dramatically influence the linear vs. nonlinear comparison. 
See also the Supplementary Information. 

We have realised a scenario proposed by Boixo et al^ to achieve metrological sensitivity 
beyond the Heisenberg limit 5X oc A^~^ using metrologically-relevant interactions among parti- 
cles. To generate pairwise photon-photon interactions, we use fast nonlinear optical effects in a 
cold atomic ensemble and measure the ensemble magnetisation ( ) with super-Heisenberg sen- 
sitivity SFz (X A^"^/^. To rigorously quantify the photon-photon interaction and the sensitivity, 
we calibrate against a precise, non-destructive, linear measurement of the same atomic quantit)!^, 
demonstrate quantum-noise-limited performance of the optical instrumentation, and place an up- 
per limit on systematic, i.e., non-atomic, nonlinearities at the few-percent level. The experi- 
ment demonstrates the use of inter-particle interactions as a new resource for quantum metrology. 
While possible applications to precision measurement will require detailed study, this first exper- 
iment shows that interactions can produce super-Heisenberg scaling and improved precision in a 
quantum-limited measurement. 
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Methods 

Linear & nonlinear probe light. The probe beam is aligned to the axis of the trap with a waist 
of 20 |Lim, chosen to match the radial dimension of the cloud. In the linear probing regime we use 
a train of forty 1 |lis pulses, pulse period 10 |lis, each containing 3 x 10^ photons detuned +1.5 GHz 
from the (F — 1 ^ — 0) transition. The maximum intensity is 0.1 Wcm~^. The signals are 
summed and can be considered a single, modulated pulse. 

The nonlinear probe consists of a single Gaussian-shaped pulse with a FWHM of 54 ns. 
The maximum intensity of the nonlinear probe is 7 Wcm~^ for a pulse with 10 photons. Theory 
predicts a^^^ = at a detuning A = 27r(462 MHz) in free space. This is modified by trap-induced 
light shifts, and we use the empirical value Aq = 27r(468.5 MHz), which gives zero rotation at low 
probe intensity. 

Instrumental noise. The instrumental noise is quantified by measuring yai^Sy) vs. input photon 
number (= A^^ or A^nl), in the absence of atoms, to find contributions from electronic noise 
y(ei) ]\fO^ shot-noise = A^^, and technical noise oc A^^, as described in the Supplementary 
Information. We find V^^^"^ and V^^l^ are 3 x 10^ and 4 x 10^ per pulse, respectively, while the 
technical noise is negligible. The instrumentation is thus shot-noise-limited over the full range 
of A^ used in the experiment. The intrinsic rotation uncertainty of the nonlinear probe (50nl is 
calculated from the measured A(/)nl as (50nl)^ = (A(/)nl)^ — ^nl^- The correction is at most 5%. 

Instrumental linearity. The linearity of the experimental system and analysis is verified using a 
wave-plate in place of the atoms to produce a linear rotation equal to the largest observed nonlinear 
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rotation. Over the full range of photon numbers used in the experiment, the detected rotation angle 
is constant to within 5%, and SQL scaling is observed. 
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Figure 1: Atom-light interface a) Experimental schematic: an ensemble of 7 x 10^ ^^Rb atoms, 
held in an optical dipole trap, is prepared in|F=l,mi7 = l)by optical pumping (OP). Linear 
(Pi), nonlinear (Pnl), and a second linear (P2) Faraday rotation probe pulses measure the atomic 
magnetisation, detected by a shot-noise-limited polarimeter (PM). The atom number is measured 
by quantitative absorption imaging (AI). b) Spectral positions of the pumping, probing, and imag- 
ing light on the D2 transition. 
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Figure 2: Calibration of nonlinear Faraday rotation, a) Ratio of nonlinear rotation to 
linear rotation v^". nonlinear probe photon number A^nl- The data points and error bars indicate 
best fit and standard errors from a linear regression (/)nl = + const, for a given A^nl- The 
red curve is a fit with equation ([2]), showing the expected nonlinear behaviour oc A^nl, with 
some saturation for large A^nl- b) & c) 0nl correlation plots for two values of A^nl- The atom 
number A^^ is varied to produce a range of cpi, and (/)nl. Green squares: no atoms A^^ = 0, red 
circles: 1.5 x 10^ < A^^ < 3.5 x 10^, blue triangles A^^ 7 x 10^. The blue triangles are shown 
as a check on detector saturation, and are not included in the analysis. 
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Figure 3: Super-Heisenberg scaling. Fractional sensitivity 5Fz^^^ / { F^) of the nonlinear probe 
versus number of interacting photons A^nl- Blue circles indicate the measured sensitivity, curves 
show results of numerical modelling, and the black lines indicate SQL, HL, and SH scaling for 
reference. Scaling surpassing the Heisenberg limit oc N^l is observed over two orders of mag- 
nitude. The measured damage 77 to the magnetisation, shown as green diamonds, confirms the 
non-destructive nature of the measurement. Error bars for standard errors would be smaller than 
the symbols and are not shown. 
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1 Supplementary Discussion 

Interaction Hamiltonian. The atom-light interaction is described using collective continuous 
variables and degenerate perturbation theory in reference [13] of the main text. We repeat essential 
results: 

The electric dipole interaction /iint = — E • d, taken in second order perturbation theory, gives 
rise to an effective (single-atom) Hamiltonian of the form 



torial and tensorial components of the interaction respectively, and the atomic collective vari- 
able is J = X^J^^^ where the superscript (z) indicates the i'th atom and = ^/J — fy^ /2, 
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plus terms in S'o which do not alter the optical polarisation. Here a^^^a^'^^ describe the vec- 
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Supp. Fig. 1: Spectra of the terms of 
the effective Hamiltonian. Black curve and 
right axis correspond to a'^^^; Green curve 
and left axis to P^'^h The curves are dashed 
when the coefficient assumes negative value. 
Detuning [MHz] is relative to the F = 1 ^ 
= of ^^Rb D2 transition. The detun- 
ings of the two probing regimes used in the 
experiment are indicated. 



jy = yfxfy + fyfxj /2, jz = /z/2 and jo = fl/2. For our case of the F=l ground state, the j 
operators, defined starting from the angular momentum operators f , represent a pseudo-spin 1/2 
system involving the states m/ = ±1. In this representation, the measurement of F^, described in 
the main text, is equally to a measurement of 2 = F^. 



For F = 1 atoms, the fourth-order contribution (again ignoring terms depending only on S^) 



is: 



m q2 



?(0) c2 



(S2) 



Forourinputstate, consisting of vertically-polarized photons, i.e., ( {Sx,Sy,Sz) ) = (A^/2,0,0), 
we can drop all but the a^^^ and (5^^^ terms, because 1) terms in Sx and S^Sx, leave the initial state 
unchanged, 2) terms in Sy and S^Sy commute with the measured variable, giving no measurable 
signal and 3) the terms in Si make a contributions smaller than the /^^^^ term by a factor ^ aS^/aSq- 



The coefficients a^^^ and (5^^^ depend strongly on the probe frequency due to the excited 
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state hyperfine structure. For the D2 line of ^"Rb, from the F = 1 ground state, they are shown 
graphically in Supp. Fig. 1. 

We note also that is sensitive to more spin degrees of freedom than is . The popu- 
lation of the state |F = 1, m = 0), i.e., A^^ — Jq? appears in proportional to and produces 
polarization self -rotation. In contrast, has no dependence on this population, which cannot be 
detected by any linear measurement. 

Atomic State Preparation. The atomic ensemble contains up to 7 x 10^ ^''Rb atoms held in an 
optical dipole trap formed by a weakly-focused (52 |Lim) beam of a Yb: YAG laser at 1030 nm with 
6 W of optical power. The trap is loaded from a conventional magneto-optical trap (MOT) and 
cooled to 25 |liK with sub-Doppler cooling. The system has demonstrated high effective optical 
depth^^ (rfo > 50) and sub-projection-noise sensitivity^ of ^ 500 spins with the linear probe. 

The ensemble is polarised ( ) = Na by optical pumping with circularly polarised light 
resonant with the F = 1 ^ = 1 transition, sent along the longitudinal axis of the trap. Repump 
light resonant with F = 2 ^ = 2 is simultaneously applied via the 6 directions of the MOT 
beams to prevent accumulation in the F = 2 hyperfine level. A small bias magnetic field of 
100 mG is applied along the axis to preserve F^. 

Before each polarisation step, the state of the ensemble is reset to a fully-mixed state by 
repeated pumping from F = 1 to F = 2 and back, using resonant lasers from the MOT beams as 
described in Koschorreck et al}. During the reset process, about 10% of the atoms escape from 
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Supp. Fig. 2: Shot noise limited detection. 

Detected variance in Sy versus input photon 
number for short, intense pulses of the non- 
linear probe. Points are measured variances, 
curve is a fit to var(S'^) = Ve\-\- N^l, where 
Ve\ is the electronic noise. Negligible tech- 
nical noise is seen. The system is shot-noise 
limited above A^nl ^ 4 x 10^, i.e., for the 
full range used in the experiment. 



the trap, allowing measurement with different A^^ during a single loading cycle. 



Shot noise limited detection. Before the ensemble, a beamsplitter and calibrated fast photodi- 
ode are used to detect the input pulse energy Sq = Sx- After the ensemble, pulses are analyzed 
in the ±45° basis with an ultra-low-noise balanced photo-detector^^, giving a direct measure of 
Sy. Both signals are recorded on a digital storage oscilloscope, and rotation angles calculated as 
(/) = Sy/ (SxVThTv) where Tny are the measured transmission coefficients for the system optics 
(vacuum cell, lenses and dichroic mirrors to separate the dipole trap beam). Supp. Fig. [2] shows 
the noise v^*. power curve for generation and detection of nonlinear probe pulses, indicating an 
electronic noise contribution to var(S'^) of 4 x 10^ per pulse. This electronic noise is subtracted 
when calculating 5(/)nl in Fig. 3 of the main text. 



System linearity. To check for systematic errors, both in the apparatus and in the analysis, we 
repeat the experiment under identical conditions but with no atoms present. We mimic the Faraday 
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Supp. Fig. 3: System linearity. Sensitiv- 
ity results obtained as in Fig. 3 of the main 
text, but with a waveplate rotation in place 
of the atomic Faraday rotation. As expected 
the sensitivity shows SQL scaling, provid- 
ing a direct verification of the linearity of 
the equipment and method of analysis. Er- 
ror bars plotted are the standard errors of the 

10^ 10^ 10^ 

Photon number: A/nl measured rotation signal. 

rotation signal by rotating the wave-plate used to balance the polarimeter to give a signal equal 
to the largest signal seen with atoms. The measured rotation is independent of A^nl, and gives 
shot-noise scaling of the sensitivity, plotted in Supp. Fig. ([3]) over the range of A^nl used in the 
experiment. 

Modelling. We model the nonlinear rotation by integrating the Maxwell-Bloch equations in three 
spatial dimensions x = {x,y,z) plus time t. This semiclassical model describes the average 
rotation ( ) = ( 5^^^^*^ ) / ( Sx^^ ), which remains <C 1, while the quantum noise is given by 

5(t) = SSt'^/ { St^ ) ^ SSt^/N = l/y]V. 

In retarded coordinates ( = z and r = t — z/c, the field envelope f (x, r) and atomic state 




21 



/9(x, r) obey the coupled equations 

V£ = —V (S3) 

drP^'-[p,H{S)]+C{p) (S4) 

where V = &1 + dy + 2ikd(^ is the differential operator of the paraxial wave equation (PWE), k is 
the wave-number, C is the Liouvillian describing relaxation and the polarization envelope P(x, r) 
is 

P = nTr[pdJ = np, (S5) 

where n is the local atomic number density and is the dipole operator describing downward 
transitions. For the atom distribution, we take a Gaussian with FWHM 2aT^An2 and 2(7^ V^ln2 in 
the transverse and longitudinal directions, respectively: 

n(x) = Na{7v^^^(7l(7^)-^ exp[-rV4] exp[-zVai] (S6) 

where = + y^. 

We solve to first order in A^^ as follows. We identify a solution to the zero-atom equation 
P£ = as the input field Specifically, we take = SoevT{T)M{x.) where ey is the unit 
vector in the V direction, 

r(r) = TT-^/V-^/^ exp[-tV2T2] (S7) 
is a Gaussian pulse with FWHM 2rVln2, and 



M(x) = ^/--|^exp[-rVi^^ (S8) 
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Supp. Fig. 4: Population dynamics un- 
der nonlinear probing. Results of simula- 
tions for a Gaussian pulse with a FWHM of 
54ns, 5.7 million photons, a peak intensity of 
4W/cm^ and A = 462MHz, with an initially- 
polarized density matrix. Blue dashed line, 
pulse intensity (normalized). Symbols show 
total population for groups of states: Green 
circles, |1, 1) state; red squares, |1, — 1) and 
1 1,0) states (100 X magnified); purple stars, 
excited states (5 x magnified); turqouise di- 
amonds, F = 2 states (5 x magnified). 



where w'^{z) = wl{l + z^/z^), zr = nwl/X, and ip is the wave-front phase. This describes a 
Gaussian beam with effective area = / d^dy |M(x, y, 0)p/|M(O)p = tiwI/2. 



We find numerically the solution to Eq. (S4) with 8 = 8^^^^ as the lowest-order atomic 
response p^^\ A representative case is shown in Supp. Fig. |4| The evolving atomic state generates 
the field 



^0 



(S9) 



where G is the Green function for the PWE. 



The detected signal is { Sy) = {hcuZo) ^ J (fxdr Sy^"^ S^"^ + c.c, where the spatial integral 
is taken over the surface of the detector and subscripts H, V indicate polarization components. It 
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can be shown, e.g., using Green function techniques'^, that 

/ d^^dr^f )*4^) = / d^^dr^:^ (SIO) 
J 2i6:o J 

and thus 

{Sy) = {hwZ^y^^ [ (fxn [ drS^^'^'p^i^ +C.C. (Sll) 
2^6:0 J J 

while {S^) = {hwZoy^ Jd^xdr l^'^p = (ncjZo)"^|^op. 

Independently determined values for the model parameters gl.cjt.w^ and r are used, leaving 
only A^^ as a free parameter, found by fitting to the data. We note that A^^ determines the vertical 
position of the curve in Figure 3, and has no effect on the sensitivity scaling. In this sense, the 
model confirms the scaling behaviour with no adjustable parameters. 

Simulations indicate that loss of polarization in F = 1, and thus rotation signal, is mostly 
due to spontaneous decay into the F=2 ground level, as seen in Supp. Fig.|4} 

Sensitivity in time- and number-limited scenarios When time is limiting, the relevant sensitivity 
is 5F^r^/', where SF^ = { F^) S(j)/(j) as in Equation 1, and the measurement duration r is tl = 
40 |Lis or tnl = 54 ns for the linear or nonlinear measurement, respectively. The sensitivity can 
be calculated from the measured values ^(A^,) = 3.3 x 10"^, and B{Ao) = 3.8 x 10~^^, using 
50 = iV-V72, 0^/ ( ) = A{Al)/2 and 0^^/ ( ) = B{Ao)N^j2, We find 5Fi^V^/' = 
1.9 X lO^Hz-^/^A^L and 5Fi^^V^{f = 6.1 x 10^^}lz-^/^N~l^\ Given an equal number of 
photons A^L = ^NL = the nonlinear technique surpasses the linear at = 3.2 x 10^, well 
within the super-Heisenberg portion of the curve in Figure 3. In contrast, when time is not a 
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limited resource, the sensitivity-per-measurement is = 3 x 10^ N^^^^, and = 2.6 x 

lO^^A^i^L^^- Extrapolating, the nonlinear technique would surpass the linear at A^nl = 8.7 x 10'', 
which is however outside the super-Heisenberg portion of the curve in Figure 3. 
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